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^ ■ Abstract 

Q"^ . The recently proposed infrared modification of gravity through the introduction of a ghost scalar 

\ field results in a number of interesting cosmological and phenomenological implications. In this 

\ paper, we derive the exact cosmological solutions for a number of scenarios where at late stages, the 

ghost behaves like dark matter, or dark energy. The full solutions give valuable information about 
the non-linear regime beyond the asymptotic first order analysis presented in the literature. The 
generic feature is that these ghost cosmologies give rise to smooth transitions between radiation 
r~| ! dominated phases (or more general power-law expansions) at early epochs and ghost dark matter 

resp. ghost dark energy dominated late epochs. The current age of our universe places us right 
■ at the non-linear transition phase. By studying the evolution backwards in time, we find that the 

^ • dominance of the ghost over ordinary baryonic matter and radiative contributions persists back 

^ to the earliest times such that the Friedmann-Robertson- Walker geometry is dictated to a good 

approximation by the ghost alone. We also find that the Jeans instability occurs in the ghost dark 
energy scenario at late times, while it is absent in the ghost dark matter scenario. 
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I. INTRODUCTION 



An interesting proposal to explain dark matter or dark energy by an infrared modification 
of standard gravity was made recently in ^ (there is a rather long history of earlier attempts 
at infrared modifications of gravity, see e.g. j^)- The key idea is that gravity at the present 
epoch might be in a Higgs phase. This is achieved by the introduction of an additional 
ghost scalar field, 0, with the wrong-sign kinetic term. If the scalar has a vanishing vacuum 
expectation value (vev), one would obviously obtain an unstable vacuum since the kinetic 
energy would be unbounded from below. However, a stable vacuum can emerge if one allows 
for a more general kinetic term. The ghost scalar might then condense in a vacuum with 
non-trivial constant velocity in field space 



Notice that this vacuum breaks Lorentz invariance spontaneously by breaking time trans- 
lation invariance. The ghost's kinetic term is described by a kinetic function P{X) where 



is (up to a factor two) the standard kinetic term. The interesting observation is then that, 
once the ghost has condensed, it assumes the equation of state of either matter or vacuum 
energy. Consequently it represents an interesting candidate for both dark matter and dark 
energy 0. 

To study fully the cosmological consequences of the ghost condensation proposal one 
first needs to derive the corresponding cosmological solutions. Though asymptotic analysis 
at late stages, close to the condensation point, does already provide important insights, it 
cannot address the full history of the ghost's evolution. In particular one might ask what 
happens to the ghost's mass density when we adjust it, according to the proposal of jH, at the 
current epoch with the observed dark matter/energy density. Will the ghost still dominate 
at earlier times or become negligible as compared to ordinary matter resp. radiation? How 
is the scale-factor of the universe affected by the ghost at earlier times? Moreover, with the 
current age of the universe, as given by its Hubble time Hq^ = (10~^^eV)~^, has the ghost 
already had sufficient time to condense, or are we still far away from this critical point. 
Moreover, one would like to know when the asymptotic behavior sets in and can be trusted 
as an approximate solution. 

Our goal in this paper is to provide these ghost cosmological solutions for a number of 
cases, as a starting point for more detailed cosmological and astrophysical investigations of 
the ghost condensation proposal. We will see that a ghost cosmology leads generically to a 
universe which evolves from a fractional power-law expansion (for the simplest model with 
quadratic P{X) one finds a t^/^ radiative phase; for models based on cubic, quartic or higher 
power P{X) more exotic fractional power time dependences arise) to a late time matter or 
dark energy dominated phase. A ghost cosmology defined through a particular P{X) will 
therefore smoothly interpolate between two types of standard cosmologies. This in itself 
seems to be of technical interest, given that normally when going from one epoch to the next, 
one simply glues standard cosmologies together in a continuous but non-smooth fashion. 
Furthermore, we will see that at the current epoch, if the ghost indeed accounts for a sizeable 
contribution to dark matter /energy, neither of the asymptotic late or early time limiting 




(1) 







X = -d^<j)d^<P 



(2) 
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cosmological solutions applies. Instead, we find ourselves right in the highly non-linear 
transitory stage between early and late time asymptotics. Hence the ghost condensation 
process, which is shown to endure for an infinite amount of time, is far from being finished. 
Via the ghost cosmological solutions the ghost's mass density will depend in a non-trivial way 
on cosmic time. In particular, it will turn out that if the current ghost's mass density is of the 
size of the presently observed dark matter /energy density, then the ghost will always have 
dominated other contributions like ordinary matter or radiation in the past and will continue 
to do so in the future. Thus as a good approximation one can neglect ordinary matter and 
radiation at other epochs as well and get a sufficiently accurate cosmological description 
by just considering the ghost-gravity system alone. Our thus derived ghost cosmologies are 
therefore relevant not only at the current epoch but also at earlier and later epochs. One 
interesting aspect of these ghost cosmological solutions concerns the coincidence problem. 
In general, since the ghost's contribution to the dark energy density does not stay constant 
in time one would expect not only the dark matter density but now also the ghost's dark 
energy density to increase towards earlier times. This should at least lead to some alleviation 
of the coincidence problem and would be interesting to study in future work. 

Finally, it would be nice to embed the ghost condensation proposal into M-theory. The 
ghost, which appears only derivatively coupled, satisfies a shift symmetry and could therefore 
easily be identified as an axion. The hope would be that M-theory can fix the magnitude of 
the ghost dark matter or energy which in the effective field theory framework needs to be 
adjusted to the observed values. Moreover M-theory predicts dark matter/dark energy on its 
own (see e.g. 0) and it would be interesting to see the interplay of these contributions with 
the ghost contribution. For other interesting implications of the ghost condensate proposal 
or related subjects, see the references in 

Without further a due, we will start the next section by presenting the general framework 
for cosmology with the ghost scalar. Section III will be devoted to the full cosmological 
solution related to the appearance of ghost dark matter at late epochs. Here, as we have 
alluded to already, we find a transition from a radiation dominated early epoch to a ghost 
dark matter dominated late epoch. Section IV gives the cosmological solution for P{X) = 
X" which arises in any Taylor-expansion of a more complicated kinetic function and is 
in particular relevant at early times. The outcome is a simple fractional power-law time 
dependence. Section V derives the cosmological solutions which give rise to a late time 
(ghost) dark energy accelerated expansion. The general result is that the full solution 
interpolates between an early epoch radiation dominated phase and the mentioned (ghost) 
dark energy dominated late epoch. In Section VI, we discuss the possibility of a transition 
between an early period with large vacuum energy, suitable for instance to sustain infiation, 
and a late period with a smaller amount of dark energy. Section VII includes ordinary 
baryonic matter and radiation in addition to the ghost. We study their evolution and find 
that if the ghost accounts for all of dark matter resp. dark energy observed today that the 
dominance of the ghost will exist not only today but persists back to much earlier times. 
For the evolution of the scale-factor of the Friedmann-Robertson- Walker (FRW) geometry 
one can therefore neglect other sources besides the ghost to a very good approximation. The 
final section VIII analyzes the potential Jeans instability of our cosmological solutions. 
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II. GHOST COSMOLOGY 



In this section, we will present the general framework describing the cosmological evolu- 
tion of a FRW spacetime in the framework of general relativity coupled to a ghost scalar. 
In later sections, we will then apply this formalism to concrete cases covering late time 
generation of ghost dark matter or dark energy. 



A. The Cosmological Equations 

Following the proposal of Q, one replaces the usual kinetic term X for the ghost scalar 
through some more complicated function P{X) leading to a ghost gravity theory defined 
by an action 

Notice that it is not necessary to add a cosmological constant term separately by hand as it 
can be absorbed into a constant added to P{X). We will elaborate more on this issue later 
in section V. The mass scale M is not determined a priori by the theory itself but should be 
regarded as a free parameter which needs to be fixed by phenomenological considerations. 
The range of M becomes thus confined to Q 

ImeV <M< lOMeV (4) 



where the upper bound can possibly even be relaxed to the GeV regime By variation of 
this action w.r.t. the metric and 0, we obtain the Einstein equation 

Gfj,u = SttGT^^^u (5) 

with ghost energy-momentum tensor 

T^,^, = M\P{X)g^, + 2P\X)d^^d,(t>) (6) 

and the equation of motion for 

d,{V^P\X)d^<P) = . (7) 

For the ghost condensate proposal, it is essential that P{X) exhibits a minimum at some 
positive value C such that X becomes timelike at this value and consequently acquires a 
non-vanishing time- dependent vev in this vacuum state. A stability analysis of fluctua- 
tions around the condensate point implies the constraints 0| 

P'{C) > , P\C) + 2CP"{C) > . (8) 

Besides these requirements, P{X) is a priori arbitrary. Though it will soon turn out that the 
Friedmann equation (for the case with flat spatial sections studied in this paper) imposes a 
further consistency constraint, it is typically the case that this additional constraint is, for 
the cases of P{X) studied in this paper, satisfied wherever (jSI) holds true. 

The ghost condensate breaks time translation invariance and the shift symmetry in 
down to a diagonal subgroup |H while leaving all other spatial symmetries like e.g. rotation 
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invariance intact. A spatially homogeneous and isotropic universe is therefore compatible 
with the idea of a ghost condensate and will be adopted throughout this paper. The corre- 
sponding geometry is described by an FRW metric where our focus is on the case with flat 
3d sections 

ds^ = -dt^ + a\t) {dr^ + r\de^ + sin^ ed(l)'')) . (9) 

In compliance with the homogeneity and isotropy assumption, the ghost scalar cannot de- 
pend on the spatial coordinates but only on time 

= m ■ (10) 

We therefore see that in the FRW background X depends only on the ghost's time derivative 

X = ^'. (11) 
Consequently the equation of motion for simplifies to 

|(a3(t)0P'(X))=O, (12) 

which is solved by 

a^{t)^P'{X) = const . (13) 

If at late times we want to have an ever increasing scale-factor a ^ oo then either or 
P'{X) have to vanish. The former choice, which would amount to redshifting the ghost 
velocity to zero, does not take us to the ghost condensate vacuum state and will not be of 
interest to us here. We will therefore adopt the second boundary condition, namely that 
asymptotically, 

X with P'(C) = as t ^ oo , (14) 

which describes the condensation of the ghost (f){t) ±\/Ct. 

This boundary condition implies that the condensation process takes place in an infinite 
amount of time. In view of the fact that the condensation point X = C for the simplest 
quadratic kinetic function P{X) = |(X — C)^ marks the border between stable and unstable 
states [Ij this is what we want. Otherwise, one would have to be concerned about 'over- 
shooting' into the unstable regime when condenses. Nevertheless, one might also wonder 
whether it might be possible for the condensation process to finish in a finite amount of 
time. If we assume this to be true then we learn from Eq. ()13|) that a{tc) = oo would have 
to hold at finite time tc indicating that the solution possesses a singularity which can be 
reached in finite time and is geodesically incomplete. We will discard such solutions and are 
hence forced to assume that the condensation process terminates at tc = oo. 

From the Einstein equation, two independent equations are obtained. Thus, it might 
seem that the problem might be overdetermined: three equations (two from the Einstein 
equation plus the ghost equation of motion) for just two unknown functions (the scale-factor 
a{t) and the ghost 0(t)). This is, however, not true. To this end, let us note that 

^"^Mt = -S'Ma'VxP') , (15) 



5 



where we have used that P = P'X. Hence, local energy conservation is equivalent to the 
equation of motion for the ghost given earlier in Eq. (|T2|l . 



V^T^^^ = ^ EOMforc/). (16) 

Since local energy conservation is already implied by the Einstein equation, the same is 
true for the ghost equation of motion. We are therefore left with just two independent 
equations coming from the Einstein equation for the two unknown functions a{t),X{t). By 
noticing that the ghost depends only on the time coordinate, we can always replace its time 
derivative by a/X and thus work exclusively with the more convenient variable X. This 
would no longer be possible if we allow for an inhomogeneous or anisotropic universe as in 
these cases, the ghost would depend on more than just one coordinate. 

Having dealt with the scalar equation of motion, let us now consider the equations for the 
metric. For the FRW background, the Einstein equation leads to two independent equations, 
the Friedmann equation 

2 

Ttl 

= —{2XP'{X) - P{X)) (17) 
3 



and the evolution equation 



H^ + 2- = -m^P{X) (18) 



where, as usual, the Hubble parameter H denotes a/ a and we have defined the mass- 
parameter m through 

m = VS^M^ = —— , (19) 
Mpi 

with Mpi the reduced Planck mass. Like the "see-saw" relation of Grand Unified Theories 
employed to explain the smallness of the neutrino masses [1], we have here a similar "see- 
saw" relation where M is the geometric mean of m and Mpi. With a range of M between 
ImeV and IGeV the parameter m lies between the extremely small energies 

4 X IQ-^^eV < m < 4 X 10~%V . (20) 

Notice that the lower bound corresponds to the current Hubble scale which is Hq ~ 10~^^ eV. 

Combining the evolution and the Friedmann equation to eliminate H from the former 
yields the acceleration equation 

^ = -^iXP'{X) + P{X)) (21) 
which is more convenient to use than the evolution equation. 

B. The General Formal Solution to the Cosmological Equations 

The Friedmann equation (for the 3d flat case with A; = considered here), due to the 
positivity of H"^, imposes another consistency constraint 

2XP'{X) > P{X) , VX(t) (22) 



6 



on any candidate ghost kinetic function P{X). We will see shortly that this constraint is 
equivalent to demanding a non- negative gravitating mass density p,^ coming from the ghost. 
It should, however, be noted that further sources of positive mass densities or a negative 3d 
curvature would contribute positively to H^, thereby relaxing the above constraint some- 
what. Let us now proceed by formally solving the Friedmann and the acceleration equation. 

First, one observes that the Friedmann equation, a first order differential equation for 
a(t), can be formally solved in terms of X through (a, = a{ti)) 



a{t) = a,e^(*'*') , f{t, U) = ±^ /* dt' ^2X{t')P'{X{t')) - P{X{t')) 

v3 Ju 



(23) 



with ti an initial reference time. Now, we can plug this solution into the Ihs of the acceleration 
equation to eliminate the scale-factor and obtain the following first order differential equation 
for X{t) 



2(P'(X))X + P = T2V3rnXP'{X)^/2XP'{X) - P{X) 
which can also be written more compactly as 
d 



— \n\{P'{X))'^X\ = T'2V3m^/2XP'{X) - P{X) . 

(JjO 



(24) 



(25) 



This is then the equation, which when supplemented with the boundary condition ()14p. 
determines X{t). Assuming that X{t) has been found by solving (j^D) we can then evaluate 
explicitly the function f{t,ti) and thus determine the scale-factor from ()23|) . Indeed, with 
help of (j2Hl) we can show for general P{X) that 



f{t,U 



i/.t'Aln|X(PW)^hln 



x,iP'{x,)Y 

X(t)(P'(X(t))^ 



1/6 



(26) 



which gives us for the scale-factor 



a{t) 



X.{P'{X,)f 



x{t){P'{x{t)y 



1/6 



(27) 



As can be easily verified, this solution for a{t) satisfies ()13p and therefore is a solution to 
the ghost equation of motion. This is in accordance with what we said earlier about the 
redundancy of the ghost equation of motion as it is already implied by the Einstein equation. 

In the following, we will solve for three classes of functions P{X) and study the 
corresponding ghost cosmologies. But before doing so, let us briefly return to the energy- 
momentum tensor In general the diagonal components of the energy-momentum tensor 
represent the mass density and pressure, T^u = {p,pgrr,pgee,pg(i><j))- By comparison to (jHl) 
one may therefore also infer a pressure and a mass density associated with the ghost field 



M^(2XP'(X) - P(X)) 



M^P(X) . 



(28) 



These quantities will turn out to be very useful in understanding various limits of our 
cosmological solutions. The constraint (|22j) implies that p^f, can never become negative. 



7 



III. GHOST DARK MATTER CASE: TRANSITION FROM RADIATION- TO 
MATTER-DOMINATED UNIVERSE WITH P{X) = i(X - Cf 

Let us now apply the formalism of the previous section to some concrete models. The 
simplest choice for P{X) letting the ghost condense at some finite positive value X = C is 
given by the parabola 

P{X) = ^-{X-Cy. (29) 

The model defined by this kinetic function will be relevant for the question of ghost dark 
matter at the present epoch as we will see. Both the stability criterion (jH| and the constraint 
()22|) require that we work in the regime where X > C > 0. Notice that in the excluded 
regime of negative X, it would be required by (fTTj) that we work with a complex ghost field 
while for positive X, it is consistent to work with a real field. 
The Friedmann and acceleration equation become 

2 

Ttl 

H^ = —{X-C) (3X + C) (30) 
6 

^ = -^{X-C){3X-C), (31) 

and it follows from the acceleration equation that any cosmological solution must describe 
a decelerated expansion for which 

a < . (32) 

where equality is attained when the ghost reaches the condensation point X — > C at late 
times. 

The differential equation ()25|) determining X{t) can be written most conveniently as 

^^^^„„XiI^(piX^ (33) 

and gets solved by the following implicit expression for X{t) 

Q{Y) + 2arctang(F) - Q, -2&rct&nQi = ±V6Cm{t - U) (34) 

where 



3y(t) + 1 x(t) 
Qm = yY(^, nt) = ^. (35) 

and Qi = Q{ti) denotes the value at initial time tj. Subsequently, we will choose the solution 
with positive sign leading to a time flow in positive direction while the ghost scalar condenses. 
To determine its value we impose the following initial boundary condition 

X{ti) > C, Y{ti) > 1 as t^t, = 0, (36) 

which fixes the integration constant with 

Q^ = ^/3. (37) 
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Notice that at late times where the boundary condition must be ()14p . it cannot be determined 
as it drops out in this hmit. This solution for X{t) maps the time interval ti = < t < oo 
faithfully into the X interval oo > X{t) > C. The ghost condensation process is therefore 
seen to take an infinite amount of time. Therefore, once the cosmological evolution starts 
within the regime where X > C, we don't have to worry that we might be driven towards the 
unstable regime where X < C in finite time. The time-dependence of the FRW scale-factor 
follows from (j27|) as 

""^^^ X(t)V6(x(t) - C)V3 Y{ty/^{Y{t) - 1)1/3 • (^^) 

Before we discuss this solution, let us first look into its early and late time behavior. 
Close to the initial time where 

t^O , Y ^oo , (39) 
we can expand the solution for Y{t) and a{t) and find 

^ Tib 

a{t) oc (Cmt)^/2 _ (41) 

Hence at early times, the ghost scalar behaves like radiation and we obtain the characteristic 
a{t) oc t^^"^ behavior of a radiation- dominated type universe. 

On the other hand at late times when the ghost has nearly reached its condensation point 
{Y approaches 1 from above) 

t^oo , Y ^1, (42) 
expansion of the solution leads to the asymptotic solution 

a{t) oc {Cmtf/^ . (44) 

Thus, as expected, the ghost field behaves at late times like non-relativistic matter with 
the characteristic a{t) oc t^/^ behavior. It could therefore account for some or all of dark 
matter at the present epoch. The time evolution of this type of ghost dark matter will 
be investigated in section VII. The full cosmological solution ^^^,^EE) therefore describes a 
smooth transition from a radiation dominated early epoch to a ghost dark matter dominated 
late epoch. 

This transition from a radiation- to a matter-dominated universe can also be understood 
by inspection of the ghost's matter density and pressure p^ at early and late times. For 
t 0,X{t) > C the expressions (j2Hl), with P{X) given in (jMO, yield 

P<k (45) 

while at late times towards the end of the condensation process, t ^ tc = oo, X C we 
obtain 

= 2M'C{X -C) + 0{{X - Cf) » p^ = 0{{X - CY) . (46) 
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FIG. 1: Plot of the ghost field X{t) (red middle curve) as a function of time for a theory based 
on P{X) = ^{X — C)^. It is clearly visible how X{t) interpolates between the early time radiation 
dominated behavior (green lower dashed curve) and the late time ghost dark matter dominated 
behavior (green upper dashed curve). The transition occurs during the relatively short time interval 
y/QCmt ~ 0.1 - 10. 



Hence, the ghost field assumes indeed the correct equation of state to behave as radiation 
at an early epoch and non-relativistic matter at a late epoch. 

We will see later in section VII that with the current age of the universe to = V-f^o 
we are led to \/6CmtQ ~ 1. If we now plot the full solutions for X{t) (Fig[T)) and a{t) 
(Fig|2)) they show besides the discussed transition from a radiation to a matter dominated 
cosmology that we are currently right in the transition period. As an explanation of dark 
matter through the ghost, this is promising as we see from ()4(i|l that towards the very end of 
the condensation process will vanish and hence its contribution to non-relativistic matter 
as well. 



IV. BEHAVIOR AT LARGE X AND THE GENERAL CASE: P{X) = X" AND 
P{X) = ix - C)" 

Exact cosmological solutions for more complicated P{X) can become quite involved. It is 
therefore important to understand the general characteristics. Let us now concentrate on the 
cosmological evolution when X S> C, i.e. at early times when the ghost is still far away from 
its condensation point. Obviously in this regime, any sufficiently smooth function P{X) can 
be Taylor-expanded in terms of powers of X. Let us therefore study first the dynamics of 
just the simple power-law 

P{X) = X"" , (47) 
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FIG. 2: Plot of the scale factor a(t) as a function of time for a ghost cosmology with -P(^) = 
^(X — C)^. The scale factor (red middle curve) evolves from a t^/^ radiation dominated behavior 
at early times (green upper dashed curve) to a t^l"^ ghost dark matter dominated behavior at late 
times (green lower dashed curve). 



which would occur in the Taylor expansion of more complicated function P{X\ In this case 
the differential equation for X(t) is easily solved through 

with initial value Xj = X(tj) at initial time tj = 0. If we assume that initially 

X, > 1 , (49) 

and take the solution with positive sign, then the result is a simple power dependence for 
X(t) as well as for a{t) (by plugging X(t) into the general (|77j)) 

-w^f^r'"i (50) 



a{t) oc . (51) 

Indeed, it is easy to understand how this time-dependence arises. Namely, for -P(X) = X" 
we find that the mass density and pressure of the ghost obey 

= (271 - l)p^ . (52) 
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Expressed in terms of the usual equation of state p = wp, this means that = -2^—[- Given 
that for a FRW universe with flat spatial sections and constant equation of state parameter 
w one gets 

a{t) oc , (53) 

we see that with the solution (|HT|) immediately follows. As a particular case, we recover 
the radiation dominated universe if n = 2. 
The more complicated kinetic function 

P{X) = (X - CY (54) 

which would, in contrast to the approximate P{X) = X", correctly describe a condensation 
process and generalize the simple parabolic P{X) of the last section gives rise to the following 
equation for Y{t) 

[{2n -1)Y -1]Y = T{2V3C^m)Y{Y - 1)"^ [{2n - 1)Y + l] ^ (55) 

which has a solution in terms of a combination of two Appell hypergeometric functions. 
Since this solution is rather involved, we will not present it here but confine ourselves 
to a discussion of the early and late time behavior. During the initial stages we have 
approximately a P{X) = X^ behavior with the power-law expansion (j51|) while at late 
stages close to the ghost's condensation point when X(t) ^ C we find a ghost equation of 
state 

P± = (X- C)/2nC (56) 

indicating a matter dominated period which once more might be used to describe dark mat- 
ter. With the generalization to arbitrary powers, we therefore recognize that transitions 
from more general power-law cosmologies a{t) oc j({2"--i)/3n. early epoch to a late epoch 

ghost dark matter dominated a{t) oc t^^^ cosmology are possible as well. 



V. DARK ENERGY CASE: TRANSITION FROM RADIATION- TO DE SITTER 
UNIVERSE WITH P{X) = \{X - Cf - D 

Let us next see how one can accommodate for a dark energy dominated de Sitter expan- 
sion phase at late times within the ghost gravity framework. Adding an explicit cosmological 
constant A to the ghost-gravity action leads to [3| 



/^From this action and also from the ensuing Einstein equation 

G^, = 87rGT^,^,-A(7^, (58) 

with T^^^i^ as given in (jHl), it is obvious that a cosmological constant term can be completely 
absorbed into P{X) simply by adding a constant 

P{C) = - A (59) 
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to the kinetic function, P{X) — > P{X) + P{C). We can therefore still work with the simpler 
action Q and use our formalism as presented in section II. Unlike the simple parabolic 
P{X) case which we studied before in relation to dark matter which had P{C) = 0, we 
would now require a non-vanishing P{C) 7^ at the condensation point X = C which is 
characterized furthermore through P'(C) = 0. This will lead to either de Sitter (dS) or 
anti-de Sitter (AdS) universes whenever such a critical point is reached 

P{C) <0 A>0 dS (60) 

P(C) > A < AdS (61) 

The simplest theory exhibiting such a late time dS expansion arises by adding to our 
previous parabolic P{X) a positive constant D 

P{X) = ^-{X-Cr~D, (62) 

such that at the critical point P{X = C) = —D. Upon condensation the ghost then behaves 
like a positive cosmological constant 

A = m^D . (63) 
The Friedmann and acceleration equations for this specific case become 

= T ~ C) (3X + C) + d) (64) 
i.4(l(A--C)(3A--C)-fl). (65) 

It is obvious from the acceleration equation that at early times where X ^ C the ensuing 
ghost cosmology undergoes a period of deceleration, a < 0, while at late times when the 
ghost condenses and we have X C this will turn into a period of acceleration, d > 0. 
Working once more in the regime where X > C will ensure that both constraints, (jH)) and 
fl^ . are satisfied such that the vacua considered are stable and the ghost mass density is 
always non-negative. This regime is also compatible with the boundary condition (|T^ at 
late times. 

To find the corresponding cosmological evolution we have to solve the differential equation 
for X{t) (j2H)). To this end, it is most convenient to work with the rescaled variable Y{t) = 
X(t)/C and to introduce the parameters 

2D 

771 = ±\/6mC , d = — - . (66) 
The differential equation then becomes 



Y{3Y - 1) = -mY{Y - 1) ^/{Y - 1){3Y + 1) + d . (67) 
Once a solution to this equation has been found, the scale-factor follows from p7jl as 

«(^) - y(t)i/a(y(t) - i)V3 " (^8) 

We will now present the solutions in different classes according to the value d which deter- 
mines the late period cosmological constant. 
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A. d>l,d^j 



For the first case d > 
energy the differential equation yields the solution 



1, (i 7^ I which comprises the case with arbitrary large vacuum 



mt 



1 



Vd 



:ln 



u±it) 



U. 



±,i 



Vd 



In 



V{t) 



V, 



where 



u±{t) 

V{t) 



Y(t) 



±{d-l-Y{t)) 
d+2{Y{t) 



(Yit) - 1) 



+ v/3r2(t) - 2Y{t) 
- 1 



1 + d 



Vd 



+ ^/3Y^{t) - 2Y{t) -1 + d^ 



(69) 

(70) 
(71) 



with initial values U±^i, Vi at time ti = and either both plus signs or both minus signs 
have to be chosen. For d = 4/3 the function U±{t) would vanish identically in the relevant 
interval 1 < Y{t) < oo. We will therefore study this case separately and concentrate here 
on the case with d ^ 4/3. Moreover, we will focus on the solution with plus sign choice, 
i.e. U+. 

The first thing to notice is that both U±{t) and V{t) have no zero in the interval 1 < 
Y(t) < oo and therefore do not change sign. Moreover, both \U^(Y(t))\ and V{Y{t)) 
decrease monotonically with Y over the interval 1 < F < oo and cover the range 



oo > ViY) >— + V3, 
y/d 



d-2 
Vd^ 



Vd 



>|f/+(>^)l>2 



1 



Vd^ 



V3 



(72) 



We therefore learn that \m\t decreases monotonically with Y over the interval 1 < F < oo. 
Hence, time increases if the universe evolves from large to small Y values. So if we let the 
universe start off at initial time = at some value Yi > 1 and condense at later time tc at 
value Y{tc) = 1, then this evolution is described by 



\m\t 



Vd^ 



■An 



+ ^In 

Vd 



Vit) 



V,. 



(73) 



By inverting this implicit solution one obtains the scale-factor via (j27j) . Notice that at the 
condensation point t goes to infinity, thus 



tr 



oo 



(74) 



and we see again that it takes an infinite amount of time for the ghost to condense which is 
why we do not have to worry about entering the unstable regime where X is smaller than 
C. 



B. d=l 

For the special case with d = ^ which we have left out in the previous analysis one finds 
the simple explicit solution 

r(t)=(l-e-^(l--)j (75) 
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with value 1^ > 1 at initial time ti. To have a solution in the stable regime F > 1 we choose 
the positive sign for m. Again, we see that the condensation process takes an infinite amount 
of time until it reaches Y{tc = oo) = 1. With ti = the ensuing scale factor becomes 



\rh\t 

a(t) (X esVs 



\rh\t 
' V3 



1 



(76) 



d=l 



For the case with d—1 the solution reads 

'W{t) + 1\/ Wi-1 



fht = 2 In 



V + 1 )\W{t)-l) 



{W{t) - Wi) , 



where 



w{t) = 



/I + 3(3r(t)2 - 2Y{t)) - 1 
Vl + 3(3F(t)2 - 2Y{t)) + 1 



1/2 



(77) 



(78) 



and initially at tj = we have W{ti) = Wi. 

Notice that the interval 1 <Y(t) < oo is mapped bijectively into the interval 1 <W{t) < 
\/3 where Y(t) = 1 is mapped to the logarithmic branch point at W{t) = 1 while Y{t) = oo 
becomes W{t) = \/3. Again the evolution of the ghost cosmology starts off at some value 

> 1 at initial time ti — and condenses towards Y{tc) ^ 1 at later times. Since Y{tc) — 1 
gives a logarithmic singularity in the solution, we see that the condensation process again 
takes an infinite amount of time, i.e. we have to set t^, = oo. Moreover, we have to choose 
the solution with positive m to have a time-flow in positive direction. If the initial condition 
is such that Yi^ 1 the solution becomes 



\m\t = 2 In 



1 



v^- 1 

W{t) V3 + 1 



( W{t) 



W{t) + V3 



(79) 



D. d<l 

Finally for the last case with d < 1 one finds 

' ' ^'W (80) 



1 / R{t) \ / Ri \\ 2 , 
\m\t — , arctan , — arctan , H ;= m 



V, 



(sending \m\t — >• — |m|t also gives a valid solution, however with opposite time arrow) where 
we have defined 

^ ' ^W^{t)-2Y{t)-l + d ^ ' 

and initial values Ri^Vi at ti = 0. It is easy to check that R{Y{t)) is a monotonically 
decreasing function along the interval 1 < y < oo covering the range 

^ ^'^<« > Tz ■ '''' 
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Therefore, the time t is monotonically decreasing with Y . Thus, with time running forwards 
we have to evolve from large to small Y values. At the condensation point, l^(tc) = 1, we 
have K(tc) = oo from which we see that once again the condensation process of the ghost 
will take an infinite amount of time, hence tc = oo. 

E. Transition from Radiation to Dark Energy Dominated Universe 
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FIG. 3: Plot of the ghost field X{t) (upper red curve) as a function of time for a ghost cosmology 
with P{X) = 1{X-Cf-D; D = fC^. One can clearly see the transition from an early radiation 
dominated universe (steep green dashed line) to a late time (dark) vacuum energy dominated 
universe (flat green dashed line). The transition occurs during ^/GCmt ~ 0.1 — 5. 

To make these cosmo logical solutions more transparent, let us look at their behavior at 
early and late times. By expanding the solutions around Y{ti) = Yi ^ 1 ai early times 
t — > = one finds for all the different cases that 

X{t) = (83) 
^ ^ V2MH ^ ^ 



a(t)oc^^tJ . (84) 

On the other hand at the late stages of the evolution upon condensation of the ghost field 
where Y{t) 1 and t ^ tc = oo one finds for all cases the general behavior 

X{t) =c(l + e'^^^A (85) 



D M 



2 



a(t) oc 3 Afp,\ (35) 
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Here we have expressed these general results in the original parameters. Hence all the 
obtained solutions describe irrespective of the value of Z) a transition from an early radiation 
dominated epoch to a dark energy dominated de Sitter cosmology at late epoch. At the same 
time this means a transition from deceleration to acceleration as stated at the beginning of 
this section. 
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FIG. 4: Plot of the scale-factor a{t) (red middle full curve) as a function of time for a ghost 
cosmology based on P{X) = ^{X — C)^ — D; D = |C^. It interpolates during the time interval 
\/6Cmt ~ 0.2 — 7 between a t^^^ behavior (radiation dominated, lower green dashed curve) at 
early times and a eyip{y^ D/3mt) behavior (vacuum energy dominated, upper green dashed curve) 
at late times. 

Again the behavior of these limiting cases can be quickly understood by inspection of the 
density and pressure associated with the ghost field and the related equation of state. From 
()28|) we obtain that at early times 

P<p , (87) 

while at late times 

P<p -P^ ■ (88) 

These are indeed the right equations of state showing that the ghost field behaves like 
radiation at an early stage and like vacuum energy at a late stage. 

Let us now look in more detail at the complete transition described by the given full 
solutions. Since all these solutions show a very similar dependence, we will pick for con- 
creteness the easiest case of d = 4/3 and choose its initial value 3> 1 which leads to a 
late time (dark) vacuum energy, M^D = ^M'^C^. The explicit transition from the radiation 
dominated to the vacuum energy dominated phase is depicted in Fig 01 for Y(t) and in Fig 01 
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for the scale- factor a{t). One observes that the transition occurs during the time interval 
y/QCmt = 0.1-5. We will see in section VII that the current age of the universe estimated 
hy to = 1/Hq leads to \/6Cmto ~ 1. The plots in FigOland FigElshow that this is right in 
the middle of the transition phase - very similar to what we have found before in section III 
in the dark matter case. 

VI. FROM LARGE TO SMALL COSMOLOGICAL CONSTANT? 

Both theories with quadratic function P{X) which we have studied in section III and 
V enjoyed a radiation dominated behavior with characteristic a{t) ~ t^/^ at early times. 
Moreover we have seen in section IV that with a higher power P{X) we can also generate 
more general but exotic power-law expansions at an early epoch. In view of the inflationary 
phase expected to have happened in our universe at very early times, it will be also interesting 
to ask whether such an inflationary phase could be described by a suitable P{X). Moreover, 
in view of the current astronomical evidence for dark energy P], one would like to join it 
smoothly after various intermediate epochs again to a late time accelerated de Sitter type 
expansion, however this time with a drastically smaller vacuum energy. Irrespective of the 
intermediate epochs which might differ from a de Sitter type expansion, let us now - by 
assuming that the ghost cosmology framework applies way back to the times of inflation 
- analyze the structure of the corresponding P{X). To describe a transition between a de 
Sitter phase with large cosmological constant and one with a small cosmological constant 
the requirements on P{X) are: 

• P{X) must possess (at least) two extrema, X = Ci and X = C2 with P'{Ci) = 
P'{C2) = at negative P{Ci) < P{C2) < 0. The latter condition guarantees a positive 
vacuum energy at the extremal points (to realize a very small late-time cosmological 
constant, one would furthermore require a suitably small vacuum energy — P(C2)M^ ~ 
meV^). 

Given such a P{X) we must moreover require that 

• The cosmological evolution must run from Ci - the extremal point with larger vacuum 
energy - to C2, the point with smaller vacuum energy. This guarantees a transition 
from a large positive cosmological constant Ai = —P{Ci)M'^/Mpi at early times to a 
small cosmological constant A2 = —P{C2)M'^ /Mpi at late times. 

Two explicit classes of functions satisfying the first requirements on P{X) are the cubic 
(with positive parameters b,D > 0) 

P{X) = -{X - + b{X -Cf -D (89) 

for which the critical points are 

Ci = C P{Ci) = -D (90) 

02 = 0+^^1 p(C2) = -D + ^b' (91) 

and the Mexican hat kinetic function (with positive parameters g,D > 0) 

P{X) = -{X - Cy + -^(X -C)^-D (92) 
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whose two relevant critical points (the left minimum and the maximum in between the two 
minima) are 

Ci = C-g ^ p^Ci) = -D-^g' (93) 

C2 = C p(C2) = -D. (94) 

Both critical points are such that along [Ci,C2] the slope P'{X) is positive and also the 
remaining stability (jHJ and positivity constraints (j^^ can be satisfied along this interval. 
A more general P{X) describing various intermediate epochs might however be much more 
complicated. 

We will, however, argue that whatever the P{X) satisfying the above first requirement 
might be, no transition from large to small cosmological constant is possible within the 
present framework (based on a 3d fiat metric Ansatz 0). In order to see this, let us start 
from the differential equation for X and observe that In |(P')^X| becomes — oo at both 
Ci as well as C2. Hence, since the time evolution must bring us from Ci to C2, the time- 
derivative of ln|(P')^X| has to start off at Ci with a large positive value and end at C2 
with a large negative value. Therefore, the Ihs of has to change sign in the interval 
[Ci,C2]. Let us now distinguish two cases. First, if P'{X) > throughout the interval 
[Ci, C2] then \/2XP' — P > on the interval and thus the rhs of cannot change sign as 
would be required by the Ihs. Thus no solution of exists for the whole interval [Ci, C2]. 
Second, in case that P'{X) < at some point in the interval, (jH} would be violated and the 
cosmological evolution from Ci to C2 would suffer instabilities against fluctuations along its 
way. We can therefore conclude that in both cases no consistent solution exists which would 
describe an evolution from Ci to C2 and thus a transition from large to small cosmological 
constant. Notice that this argument is valid for general P{X) and thus cannot be avoided 
by introducing a bunch of intermediate epochs - as long as the whole evolution from Ci to 
C2 is dominated by the ghost. 

VII. EVOLUTION OF DARK ENERGY/MATTER IN THE COMBINED GHOST 
- RADIATION - MATTER COSMOS 

Let us in this section study the cosmological evolution in the presence of the ghost as well 
as additional ordinary matter and radiation described by perfect fluids. With n additional 
perfect fluids the Einstein equation governing the behavior of the FRW metric becomes 

n 

G^, = SnG{T^,^, + J2 , (95) 

i=l 

with T^,^y the ghost energy-momentum tensor (jH)) and 

Ti,^,u = {pi + Pi)u^Uy + piQf,^ (96) 

the energy-momentum tensor for the ith perfect fluid. A fluid which is comoving in the 
FRW rest frame possesses a velocity ■u'^ = (1, 0, 0, 0) such that its energy-momentum tensor 
becomes simply Ti^^,^ = {pi,pigrr,Pigee,Pig(i><j>)- 

The Einstein equation implies the local energy-conservation of the combined ghost-fluid 
energy-momentum tensor. Motivated by the tiny currently observed amount of 4% ordinary 
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baryonic matter as opposed to 70% dark matter [o!], we will assume in the sequel that the 
back-reaction of the perfect fluids on the geometry is negligible as compared to the back- 
reaction of the ghost which we will hold responsible in this section for all of the current dark 
matter resp. energy. That this assumption, once true at the current epoch, remains valid at 
earlier epochs will be shown explicitly below fv?\. It is thus reasonable that the ghost and 
the fluids satisfy separate energy conservation laws 

y^T^,^, = , V^T,,^, = . (97) 

As well-known, the fluid energy conservation law in an FRW background is equivalent to 

p, + 3H{pi+pi) = . (98) 

Together with an equation of state 

Pi = WiPi , Wi = const (99) 

for the ith fluid one obtains for the mass density the solution 

"'W - f "» , (100) 



pi,o ^ait) 

with pifi, ao the present values. Similarly the ghost energy conservation governs the ghost's 
dynamics as it is equivalent to the ghost's equation of motion (fT6|l and is solved by fT3|) . Of 
course the local energy-conservation of the total mass density is just a consequence of the 
Einstein equation and therefore not independent. 

The Einstein equation in the FRW background provides us with two independent equa- 
tions, the Friedmann equation 



and the evolution equation 



If we eliminate in this equation with help of the Friedmann equation we arrive at the 
acceleration equation 



a 



a 6Mpi 



^ (p<^ + Sp^ + 5^(p, + 3p,)) . (103) 



The formal solution to the Friedmann equation is given by 

a{t) = aoe^(*'*°) , g{t,to) = ±^i- f dt'{p^ + J2p^f' ' (^O^) 

V SMpi J to j 

Plugging this into the acceleration equation to eliminate the scale-factor leads to the equation 

V ^ A 
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which determines the evolution of the densities once the equation of states (jSHI) are used 
together with (jlOOj) and ()104j) . The resulting differential equation is no longer straightfor- 
ward to solve. Fortunately, this is not needed since a potential explanation of dark matter 
or dark energy at the current epoch through the ghost would imply that today the ghost 
dominates the fluid components by far and we will shortly see that this feature also persists 
if we evolve the ghost-radiation-matter system back to much earlier and also later times. 

Nevertheless, before proceeding with the full combined ghost-radiation-matter system 
let us look at two important limiting cases. The first is the limit of early times for the 
theory given by P{X) = ^{X — C)^ — D. At early times we know that the baryonic matter 
contribution is negligible as compared to radiation. The latter has an equation of state 
Pr = "^Pr- Moreover at these early times we have found that for this P{X) the ghost also 
behaves like radiation with = 3p^. Therefore in this regime (jlOSj) can be easily solved by 
(choosing the relevant minus sign) 

p^ + p^ = J (106) 

leading via ()l()4j) to 

t ^ : a{t) oc t^/^ . (107) 

The second interesting limit is for late times for the theory defined by the kinetic function 
P(X) = \iX — C)^. This kinetic function did not lead to dark energy at late times but 
instead to a non-relativistic dark matter behavior with ^ p^- Furthermore we know 
that at late times ordinary matter dominates over radiation such that we have here for the 
perfect fluid an equation of state pB ^ Vb as well. Also in the late time limit ()105|) is easy 
to solve through (again with the relevant minus-sign) 

t^t,: p^ + pj, = —ZL (108) 

which via ()104|) leads to the corresponding late time scale-factor behavior 

t-^t^ : a{t) oc t^/^ . (109) 



A. Time-Dependence of Ghost Dark Matter, Radiation and Matter 

Let us now study the time evolution of ghost dark matter, baryonic matter and radiation. 
We will therefore take the theory with P{X) = ^{X — C)^ in which the ghost led to non- 
relativistic matter at a late epoch. We know from observation, that nowadays baryonic 
visible matter contributes 4% and radiation only 0.008% to the total mass density of the 
universe as compared to 26% coming from dark matter We will here equate the latter 
with the ghost dark matter. It will therefore be a good approximation at the current epoch 
to neglect the backreaction of ordinary baryonic matter and radiation against that of the 
ghost. In this ghost dominated background we will now evolve baryonic matter and radiation 
according to the simple power-law behavior (jlOOp with equation of state parameters wb = 
and wr = I and scale-factor a(t) given by ()38|) . We will shortly see that the dominance 
of the ghost over baryonic matter and radiation becomes even more pronounced when we 
look back in time. Also at late times the dominance of the ghost will be preserved. This 
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FIG. 5: Time dependence of the mass densities P(f),PB,PR for the ghost scalar (which is assumed 
to accomit for the present dark matter), baryonic matter and radiation in a ghost cosmology based 
on P{X) = ^{X — C)^. If the ghost accounts for the dark matter today then it dominates ordinary 
matter and radiation not only today but even more so in the past as well. 



will then establish the relevance of our cosmological solution (jMj) . (jHHj) which took only the 
backreaction of the ghost into account. 

We still have to fix the relative normalizations of the densities which we will do at the 
current epoch. To this end, let us estimate the age of the universe by the current Hubble 
parameter, = l/-f^o — l/(10~^^eV). To obtain a current dark matter density of order meV^ 
we have to adjust the current ghost mass density = 2C'^M'^{Yq — 1) at this order (here 
Yito) = Yq). It is therefore natural (avoiding a finetuning of Yq) to set CM^ ~ meV^. With 
this information we arrive at a value for mto = ^/QCmtQ which is close to 1. Consequently 
Fo = 2.32 and the time-dependence of the scale-factor (|HH)) becomes 

^ ^ Y{ty/%Y{t) - 

a{t) 1.26 ■ ^ ^ 

Combining this information with the present measured densities j^, (the present ra- 
diative contribution p/jo arises from the contribution of three relativistic neutrinos plus 
the cosmic microwave background photon s lllll and the assumption of a Hubble parameter 
Ho = lOO/ikmsec-^Mpc"^ with h = 0.71 jlj) 

P<t>fi ■■ Pb,o ■■ Pr,o = 26 : 4 : 8 X 10"^ (111) 
allows us to fix the relative normalization. Thus, the properly normalized mass densities 
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become 



(r(t)-i)(3r(t) + i) 



(112) 



P0,O 
PB 



10.27 



26 \a{t)) 



(113) 



8 X 10-3 / ao \4 



(114) 



26 \a(t)) 



with p^^o = Pdm,o the current dark matter density. Together with piO|) this gives all densities 
as functions of Y{t) and thus via the inverted solution ()34|) as functions of cosmic time t. 
The time- dependence of P(t,,PB,PR is plotted in Fig|Sl We see once more that the ghost 
behaves like radiation at early times and assumes a non-relativistic matter-like behavior at 
late times. The most important aspect of Fig|3 however, is that the ghost's domination 
over all other contributions persists to earlier epochs and becomes even more pronounced at 
earlier times. At later epochs the dominance of the ghost will be preserved as well. Our 
approximation to choose (j34j) as the background at all times, not only at the present epoch, 
proves therefore to be justified and establishes that our simple cosmological solution 
()38|) is actually the relevant one also in the presence of ordinary matter or radiation. 

To obtain the normalized mass densities fl^ = p^/ Pc and Vti = Pi/pc, obeying + 
Qb + = 1, one would have to divide further by the critical density pc = 3MpiH'^. 
Since we use the approximation that the ghost alone determines the geometry, i.e. that 
Pc = Pfj} + Pb + Pr — Pel,, it follows that pc resp. H (given by (ISHl)) are known functions of 
Y{t) as well. Thus the solution (j34j) fully specifies the evolution of all fi^, In particular 
in this approximation, we will have fi,^ ~ 1 at any time. 

B. Time-Dependence of Ghost Dark Energy, Radiation and Matter 

Let us next use the ghost field to account not for dark matter but instead for the full 
amount of dark energy currently observed. To this end, we have to take the theory with 
ghost kinetic function P{X) = |(X — C)^ — D leading to a late-time vacuum energy M'^D 
when the ghost condenses. Again, from observation, radiation and baryonic matter densities 
are negligible at the current epoch as compared to the 70% contribution to the mass density 
coming from dark energy which we equate with the current ghost mass density p^^. We are 
therefore entitled to neglect all but the backreaction of the ghost on the geometry which 
means that we have at present an FRW geometry with scale factor given by one of the 
solutions in section V. Since all of them behave similarly let us take the explicit solution 
(j76p for D = with initial condition Yi ^ 1. If the ghost field accounts for all of 
the present dark energy, which we will assume in this subsection, we have to require that 
M^D ~ M^C2 is of order meV^ 

Let us be more precise about the constants appearing in the scale-factor and the mass 
densities. To this end we will again estimate the current age of the universe by to = l/-f^o — 
1/ (lO-^^eV) which brings mto = \/6CmtQ close to 1. With knowledge of tg we can then fix 
the overall constant in the scale-factor a{t). Via (fTBj) this leads to 




(115) 
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FIG. 6: Time dependence of the mass densities p^,Pb,Pr for the ghost field (which is identified 
here with the dark energy), baryonic matter and radiation. Close to the present age of the universe, 
where VGCmto ~ 1, we find that the ghost field switches from a radiative phase to a dark energy 
phase. 



Moreover from ()75|) we obtain that Yq = 2.28. Next, we fix the relative normahzations of the 
mass densities through the current observed mass densities for dark energy (which in this 
subsection will be identified with the ghost mass density), baryonic matter and radiation 

P0,o : Pb,o ■■ PRfi = 70 : 4 : 8 X 10"=^ . (116) 

This gives us the properly normalized mass densities (rf = 4/3) 

^ {(Yit)-l){3Yit) + l) + d) ^ {3Y{t) - 1)' 
P^fi 11.37 34.11 ^ ' 

-^x (118) 
70 \a{t)J ^ ' 



PR ^ 8 X 10~3 ^ / gp x4 
p^,o 70 \a{t)) 



(119) 



as functions of Y{t) where p^^ = Pde,o is the currently observed dark energy density. With 
the solution for Y{t) (j75|) we finally obtain the full time dependence for these densities which 
is plotted in FiglHl 

We see from Fig El that the hierarchy among the three contributions, which exists today 
if the ghost accounts for all of dark energy, extends to earlier and later epochs likewise. 
We therefore see that the assumption to neglect the backreaction of radiation and matter 
on the FRW geometry and to take only the backreaction of the ghost into account is a valid 
assumption at all epochs. Herein lies the importance of the cosmological solutions obtained 
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in section V. Around y/QCmt^ ~ 5 x 10~^ we reach the point where the baryonic matter 
and radiation densities coincide. At t = the ghost field has aheady adopted a radiative 
behavior as it runs with the same slope as the radiation density. At even earlier times 
radiation will dominate baryonic matter as in the standard hot big bang cosmology but 
both are still dominated by far by the ghost's mass density. Notice that while both ps and 
pR have the same dependence on the scale-factor in the dark matter case (last subsection) 
and the dark energy case (this subsection) , nevertheless their dependence on time is different 
in both cases since the scale-factor a{t) has a very different time-dependence in these two 
cases. Let us mention that for the normalized densities VL^,VLb-,^r the same remarks as in 
the preceding subsection apply. In particular we have ~ 1 at all epochs. 

The plot in FiglHl reveals another interesting aspect concerning the coincidence problem. 
The coincidence problem arises in standard cosmology when the currently observed dark en- 
ergy is explained through an ordinary cosmo logical constant. The matter density associated 
with the cosmological constant stays constant during cosmic evolution. Thus, though the 
cosmological constant would be the dominant energy contribution at present, this quickly 
changes once one goes back to earlier times. Hence, the current observation that the dark 
matter, the cosmological constant and to a lesser extent also the matter contribution are 
of the same order of magnitude present the coincidence problem. Notice that this kind of 
coincidence problem is much more harmless in the ghost cosmology framework if we attribute 
dark energy completely to the ghost. Though at late times the ghost acts like a cosmological 
constant its mass density grows towards earlier times, keeping the hierarchy among the mass 
densities currently observed more or less intact. Therefore, this hierarchy is not specific to 
a finetuned snapshot of today's universe. The real question which remains is to explain the 
hierarchy at some point in time (see e.g. [13 )• The time-dependence of the ghost's mass 
density should have interesting cosmological consequences for the verification or falsification 
of the ghost gravity proposal. 

C. Describing both Dark Matter and Dark Energy by the Ghost 

Up to now we have considered the dark matter and dark energy cases separately by 
assuming that the ghost accounts for all of dark matter or all of dark energy today. This 
was motivated by the ghost's behavior towards late times. However, we have seen that 
actually at the present epoch we are right at the transition region from the early to the 
late time asymptotics. Therefore, the ghost cosmology based on P{X) = i(X — C)^ — D \s 
indeed capable of describing both dark matter and dark energy at the same time as we will 
now discuss. 

For this, let us split the ghost's mass density into a dark matter part, Pda/, and a dark 
energy part, poE 



P4> — Pdm + Pde 



(120) 



where 



PDM = \M\X-C){'iX + C) 



Pde = M^D . 



(121) 



The associated pressures are 



Pdm = ^M\X - Cf 



PDE = -M^D . 



(122) 
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FIG. 7: Time dependence of the mass densities Pde, Pdm, Pb, Pb. for the ghost field (which is 
identified here with both the dark energy and dark matter), baryonic matter and radiation. Close 
to the present age of the universe, where VGCmto = 1.674, we find that the ghost field switches 
from a radiative phase to a dark energy phase. 

Thus the dark energy component is nothing but a standard cosmological constant. Its 
evolution is therefore trivial and can be fixed once and for all by identifying p£,E with the 
current observed dark energy density pde,o- This gives a value of Y{to) = 1.614 which 
translates into mto = VQCmtQ = 1.674. Notice that we are using a different method of 
arriving at the V(to) and fhto as we have an additional quantity to fit. We find that this 
value of rhto coincides with our previous estimates which have shown that this value is of 
order one. On the other hand, the dark matter component is the same as in subsection A, 
only that the normalization is slightly different. The splitting that we have done to match 
the current dark matter and dark energy, though an important step in obtaining a realistic 
cosmology, leads us back to the coincidence problem for the constant pde- 

VIII. JEANS-INSTABILITY 

Let us finally analyze the stability of the obtained ghost cosmological solutions in section 
III and V under gravitational collapse, i.e. Jeans instability. A universe with mass density 
p is unstable under perturbations if the perturbations have a wavelength A greater than the 
critical Jeans wave length 




(123) 



We neglect here factors of order one on the rhs which depend on the detailed geometrical 
properties of the perturbed equilibrium state and the perturbations. The isothermal speed 
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FIG. 8: Plot of cosmic time t and of the free-fall time r as functions of time for the ghost cosmology 
based on P{X) = ^{X-Cf. At ah times we have t > t indicating that no Jeans instability arises 
in the corresponding cosmological solutions of section III. 



of sound is denoted by Vg- The corresponding time scale measuring the total time of the 
ensuing collapse if A > Aj is given by the free-fall time 



SttM, 



PI 



/Up 



VP 



(124) 



Thus, as long as the cosmic time t satisfies t < t, there was not sufficient time for the 
universe to have collapsed and it will be stable. To the contrary, for t > t the universe had 
enough time to collapse and will suffer from a Jeans instability. For the ghost cosmologies 



based on P{X) = ^{X — Cy — D we have to set p 



with 



M' 



(-(X - C){3X + C) + d)= ^ ((F - 1)(3F + l) + d) 



(125) 



Let us first examine the dark matter case, i.e. the theory with D = 0. The comparison 
between t ffrom IMj) and the free-fall time r is presented in Fig|H| We see that in this case 
no Jeans instability arises at any time. To understand this better let us investigate the early 
and late time limits. For the ghost mass density these are 





oo 



2CM\X - C) 



This leads to the following expressions for the free-fall time 







oo 



m{X-C) 



(126) 



(127) 
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FIG. 9: Plot of cosmic time t and of the free-fall time r as functions of time for the ghost cosmology 
based on P{X) = ^{X - Cf - D, [D = fC^ is used for the plot). At early times we find t > t 
indicating that no Jeans instability arises. However at late times we see that there is a crossover to 
a regime with t < t. Hence a Jeans instability occurs at late times for the cosmological solutions 
of section V. 



On the other hand, the cosmic time in these limits can be found from the expressions in 
section III as 



t 



t 



oo 



V2mX 
V2_ 



(128) 



Obviously in both limits, we find that t > t. Notice that in the late time limit the ghost 
condenses, X — > C. 

For the dark energy case with nonvanishing D ^ Owe see, however, from the plot in FiglHl 
that a Jeans instability arises at late times (the actual time when this happens depends on 
D but is always in the future, i.e. later than to). Let us again look at the limiting expressions 
for this case. We find for the mass densities 







oo 



(129) 



which leads to the following free-fall times at early and late times 



t 



oo 



m\/D 



(130) 
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To compare, we find for tlie cosmic time from the expressions in section V that 



V 3Z)m, 

Therefore, at early times we find t > t and consequently no Jeans-instabihty can occur. 
However, at late times we clearly have t > t due to the logarithmic divergence in t in- 
dicating a Jeans-instability at a late epoch. Notice that the precise time, t crossover , when 
the Jeans instability will set in depends on D. E.g. in the case oi D = one finds 

that \/QC rat crossover = 15.04. Since for ^ we must recover the previously analyzed 
case which had no Jeans-instability at any time, we recognize that sending D ^ Q implies 
sending the crossover time, at which r = t crossover ■, towards infinity. Thus for small enough 

the Jeans instability can be pushed far into the future. Notice that the Hubble friction 
which is present when the instability sets in, will cure it. 

Additionally, higher derivative terms will not affect the stability. At early times, the 
higher derivative fluctuations are suppressed with respect to P{X) while at late times, 
P{X) is vanishing and hence higher derivative fluctuations might play a part. However, the 
detailed fluctuation analysis done in [H shows that we do have a stable background. 



IX. CONCLUSION 

We have derived the exact cosmological solutions for a universe composed primarily of 
a ghost condensate. It was further shown that these solutions would smoothly interpolate 
between different standard cosmologies. We find, generically for quadratic kinetic functions, 
an early-time radiation-dominated universe which transitions into a de-Sitter expansion 
phase or a late-time matter-dominated scenario depending on whether the cosmological 
constant is present or not. This could also be understood in terms of the stress energy 
tensor of the ghost which gives rise to an equation of state that clearly exhibits early-time 
radiation-like and late-time dust-like behavior or de Sitter acceleration. For higher order 
kinetic functions we find more exotic power law FRW universes. An important result is that, 
given the current age of the universe in terms of the Hubble time, we are placed right at the 
transition regime. Moreover, the cosmological solutions show that the ghost condensation 
process takes an infinite amount of time. This preserves us from entering the regime of 
unstable vacua in finite cosmic time. 

As a first step towards constructing fully realistic cosmologies, we considered the evolution 
of the universe by extrapolating from present dark energy/matter, baryonic and radiation 
densities to earlier times. We find that the ghost contribution completely dominates at all 
times and hence must be the principal driving force of cosmological evolution. Once we 
adjust for the dark matter/energy densities observed today by assuming that the ghost can 
account for these, our solutions map out the entire evolution of these densities. The other 
important consideration is the Jeans instability and we have shown that for the ghost dark 
matter case no such instability arises. In contrast, the dark energy case exhibits a late-time 
Jeans instability. It occurs when the universe has already entered the de Sitter expansion 
phase which means that Hubble friction should cure the instability. 
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